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Abstract: It is shown how a recently proposed method of cloaking is 
effective over a broad range of frequencies. The method is based on three 
or more active devices. The devices, while not radiating significantly, create 
a "quiet zone" between the devices where the wave amplitude is small. 
Objects placed within this region are virtually invisible. The cloaking is 
demonstrated by simulations with a broadband incident pulse. 
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1. Introduction 

A tremendous amount of interest and excitement has been generated by recent strides towards 
making objects invisible, not by camouflage, but by manipulating the fields in such a way that 
the cloaking device and the object to be cloaked scatter very little radiation in any direction and 
do not absorb it. Here using a new method of active exterior cloaking, described in [ 1 ] for single 
frequency waves, we demonstrate how an object can be cloaked against an incoming broadband 
pulse. To our knowledge this is the first example showing cloaking of an object against an in- 
coming pulse. It uses active cloaking devices to generate anomalous localized waves which 
cancel the incident waves within the cloaking region to create a "quiet zone", within which 
objects can be hidden. Because active devices, rather than materials, are used to generate the 
anomalous localized waves, one may superimpose the results for different frequencies to ob- 
tain broadband cloaking. Our method requires one to know the form of the incoming pulse in 
advance, since the fields generated by the cloaking devices are tailored to the incoming fields. 

Dolin [2], Kerker (3), and Alu and Engheta (4) realized that certain objects could be made 
invisible by coating them with an appropriate material tailored according to the object to be 
cloaked. A breakthrough came with the work of Greenleaf et al. (5 ], for conductivity, Leonhardt 
(6), for geometric optics, and Pendry et al. [7 ] for electromagnetism who showed that materials 
could guide fields around a region, leaving a "quiet" zone in that region within which objects 
could be placed without disturbing the surrounding field. This idea was extended to acoustics 
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11), elastodynamics 1 12], and water waves (T3J, and has been confirmed experimentally 



HMD. 



A completely different type of cloaking, which we call exterior cloaking because the cloaking 
region is outside the cloaking device, was introduced by Milton, Nicorovici, McPhedran and 
collaborators (18 19 20). They showed that clusters of polarizable dipoles within a critical 



distance of a flat or cylindrical superlens |2T)|22 23] are cloaked. Anomalously localized fields 
generated by the interaction between the induced dipoles and the superlens effectively cancel 
the fields acting on the polarizable dipoles. While larger objects do not appear to be cloaked 
(24), Lai et al. | 25 ] show that an object outside a superlens can be cloaked if the appropriate 
"antiobject" is embedded in the superlens. 

Ideally cloaking should be over a broad range of frequencies. Most cloaking methods are 
narrowband and approaches to obtain broadband cloaking can have drawbacks, such as requir- 
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ing frequency independent relative dielectric constants or relative refractive indices less than 
one (26 [27 T5] 1_6 17 1 which necessitate a surrounding medium with dielectric constant suffi- 



ciently greater than one: thus such electromagnetic cloaking could work underwater or in glass, 
but not in air or space. One proposal without this drawback is the broadband interior cloaking 
scheme of Miller |28| which uses active cloaking controls rather than passive materials. Here 
we also use active cloaking devices to achieve broadband exterior cloaking. The principle is 
similar to that of active sound control (see e.g [29, 30]), with the fundamental novelty that we 
do not need a closed surface to suppress the incident field in a region while not radiating signif- 
icantly. Another type of broadband exterior cloaking, using waveguides to guide waves around 
a "quiet zone", has recently been introduced and confirmed experimentally (3TJ. 

2. Cloaking a single frequency 

For simplicity we just consider the two dimensional case, corresponding to transverse electric 
or magnetic waves, so the governing equation is the Helmholtz equation Au + k 2 u = 0. Here 
w(x,<b) is the wave field, k = 2n/X is the wavenumber and A = 27ico/co is the wavelength 
at frequency (0 and at a constant propagation speed c$. We would like to cloak a region in 
the plane from a known probing (incident) wave m/(x,G)) supported in the frequency band 
C0o + [—5/2,5/2], where the central frequency is (Oq and the bandwidth is B. 

The key to our cloaking method are devices that (a) cancel the probing wave in the region to 
be cloaked and (b) radiate very little waves away from the devices. To give a concrete example, 
let us take for the cloaked region the disk |x| < a and assume we measure the radiation emitted 
by the devices on the circle |x| = y > a. Thus the device's field m^(x, (o) must be so that (a') 
Ud « —Ui for |x| < a and (b') Ud « for |x| = y. 

The devices can be idealized by D points xi , . . . , x# with \xj ■] = 8 and a < 8 < y so that the 
devices surround the cloaked region. Because the device's field must solve Helmholtz equation 
and become small far away, we take it as a linear combination of outgoing waves emanating 
from the source points xi , . . . , x# with the form (TJ: 

D N 

U d (x,CO)=Y< Y< *m,n^n (*|X-X m |)exp[m0 m ], 

m=l n=—N 

where is the n— th Hankel function of the first kind and G m = arg(x — x m ) is the angle 
between x — x m and (1,0). We seek coefficients b m , n so that (a') holds on points of the circle 
|x| = a and (b') holds on points of the circle |x| = y. The control points are uniformly distributed 
and at most A/2 apart on each circle. The resulting linear equations are solved in the least 
squares sense with the Truncated Singular Value Decomposition in two steps. First we find 
coefficients b m , n so that (a') holds, and second we find a correction to enforce (b') while still 
satisfying (a'), (see Appendix |A| for more details). 

3. Simulations 

We demonstrate cloaking in a regime that could correspond to Transverse Electric microwaves 
in air (neglecting dispersion and attenuation), where w(x, (6) is the transverse component of 
the electric field. For the numerical experiments we took a central frequency and bandwidth 
of 2.4GHz, a propagation speed of cq = 3 x 10 8 m/s and a central wavelength of Ao = 12.5cm. 
Simulations suggest a minimum of three devices are needed to cloak independently of the 
direction of the incoming waves. In Fig.[T]we show cloaking at the central frequency COo / (2k) = 
2.4GHz of a region of radius a = 2Ao (solid white circle). Here the devices are located 8 = IOAq 
from the origin and invisibility is enforced at a distance y = 20Ao from the origin (dashed white 
circle). The incident wave is a point source originating at x^ = (— 20,0)Ao and modulated in 
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(a): x\ (in Ao) (b): x\ (in Ao) 

Fig. 1. Wave field at the central frequency (Oq when the cloaking devices are (a) inactive 
and (b) active. Only the real part of the fields is displayed, with a linear color scale going 
from —1 (dark blue) to 1 (dark red). All fields have been rescaled by |w;(x, co)\ to remove 
the geometric spreading of the point source. 
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Fig. 2. Reduction of the scattering in percent achieved by the cloaking devices over the 
bandwidth, measured in the L 2 norm on |x| = 20 Aq. The ordinates scale is logarithmic. 

frequency by a Gaussian truncated to the bandwidth g(co) = oV27texp(— o 2 (co — (Oo) 2 /2), for 
\(D — (Do\ < B/2, and otherwise. We took o = 4/coo. The scatterer is a perfectly conducting 
"kite" obstacle (32) with homogeneous Dirichlet boundary conditions, fitting inside the cloaked 
region. The scattered field is computed using the boundary integral equation method in [32]. 

With the devices inactive (Fig. [T^), the "kite" scatters the incident field and thus can be 
easily detected. When the devices are active (Fig. [I])) they create a region with very small fields 
while being nearly undetectable from far away. Since there are almost no waves in the cloaked 
region, the scattered waves from the object are greatly reduced, making the object invisible. 
Quantitatively, the disturbance of the incident wave is 1 . 1 x 10 -4 % of the field scattered without 
the devices, as measured on the dashed white circle with the L 2 norm. We carried the same 
procedure for Nf req = 101 frequencies in the bandwidth with similar results, as can be seen in 
Fig. [2] Since the bandwidth is 100% of the central frequency, broadband cloaking is possible 
with our approach. 

Several time snapshots (computed by taking the inverse Fourier transform in time of the total 
fields) appear in Fig. [3] For full animations see Fig. [3] (Movie 1 and Movie 2), which covers 
the time interval [0, T], with T « 132ns or the time it takes for the wave to travel 50.5Aq. The 
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Fig. 3. Cloaking for a circular wave pulse. Top row: devices inactive (Movie 1 ). Bottom 
row: devices active (Movie 2). The visualization window in x is as in Fig. [T] and the scale 
is linear and relative to the maximum amplitude on the plane of the incident field at each 
time. 



devices make the incoming wave disappear when it reaches the cloaked region, and then rebuild 
the wave as it exits the cloaked region. This makes the object virtually undetectable. 

Here the cloaked region is for visualization purposes deliberately small (2Ao = 25cm radius). 
However we have successfully cloaked regions up to lOAo = 1 .25m in radius as shown in Fig.|4j 
where we take three devices with |x m | = 8 = 5a and invisibility is enforced on |x| = y = 10a. 

The point-like devices could be problematic in practice because of the @( |x — x m | ~ N ) singu- 
larities near the devices. Fortunately the point devices can be replaced (with Green's identities) 
by curves where the fields have reasonable amplitudes and where we can control a single- 
and double-layer potential (32). These curves could be the circles suggested by the contours 
\ud\ = 100max| x | =a |m;(x, co) \ (in black in Fig.[Tj). The radius of such devices for other cloaked 
region radii is estimated in Fig. [5] Since the devices do not completely surround the region to 
be cloaked, exterior cloaking is possible at least for a < 10Ao. 

A theoretical study of this cloaking method is ongoing, and we speculate it generalizes to 
three dimensions and Maxwell equations. 
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A. Finding the driving coefficients for the devices 

Let b be a vector with the (2N + l)D possibly complex coefficients b m ^ n . Denote by (resp. 
pj) the N a (resp. N r ) control points on the circle |x| = a (resp. |x| = /). The coefficients b, the 
number N of terms in the expression for and the control points all depend on the frequency 
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Kf + Wrfll/Hwfll on |x| = a \\u d \\/\\ui\\ on |x| = y 




(a): a (in Aq) (b): a (in Aq) 



Fig. 4. Cloak performance in terms of the cloaked region radius a, as measured by (a) 
+ with the L 2 norm on |x| = a and (b) ||w</||/||w/|| with the L 2 norm on |x| = y. 

Dashed, solid and dash-dotted lines correspond to co/(2n) = 1 .2GHz, 2.4GHz and 3.6GHz. 

4.5r 
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Fig. 5. Estimate of the radius of circular cloaking devices relative to the radius of the 
cloaked region a. Dashed, solid and dash-dotted lines correspond to co/(2n) = 1.2GHz, 
2.4GHz and 3.6GHz, respectively. The device radius is estimated as the largest of the dis- 
tances from a device point x m to the level-set = 100max| x | =a |w/(x, co)\. If the devices 
were perfect circles, any device radius below the dotted line (at cos(k/6)8/ol = 5^/3/2) 
would indicate that the devices do not touch, i.e. they are three disjoint devices. 



CO. Using the form for in the text, construct a matrix A of size N a x (IN + \ )D and a matrix 
B of size N 7 x (IN + l)D such that wj(pj, co) = (Ab)j and wj(pj, co) = (Bb)j. We estimate the 
driving coefficients as follows. Enforce (a'): Use the Truncated Singular Value Decomposition 
(TSVD) to find coefficients bo such that Abo ~ — [w;(pj, co)], in the least squares sense. Enforce 
(b') while still satisfying (a'): Find z as the solution to the least squares problem B (bo + Zz) « 
0. Again the TSVD can be used for this step. Here Z is a matrix with columns spanning the 
nullspace of A (a byproduct of the SVD in the previous step). 

The coefficients are then b = bo + Zz. The heuristic for N in the numerical experiments is 
N = \k(S — a/2)] , where \x] is the smallest integer larger than x. The choice of cut-off singular 
values can be used to control how well one wants to satisfy (a') and (b'). We used a fixed 10 -5 
tolerance relative to the maximum amplitude of U{ on the control points. 
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